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ABSTRACT 

We observed the X-ray bright E3 galaxy NGC 1600 and nearby members of the NGC 1600 group 
with the Chandra X-ray Observatory ACIS S3 to study their X-ray properties. Unresolved emission 
dominates the observation; however, we resolved some of the emission into 71 sources, most of which 
are low-mass X-ray binaries (LMXBs) associated with NGC 1600. Twenty-one of the sources have 
Lx > 2 X lO'^^ ergs s~^ (0.3-10.0 keV; assuming they are at the distance of NGC 1600) marking 
them as ultraluminous X-ray point source (ULX) candidates; we expect that only 11 ±2 are unrelated 
foreground/background sources. NGC 1600 may have the largest number of ULX candidates in an 
early-type galaxy to date; however, cosmic variance in the number of background AGN cannot be 
ruled out. The spectra and luminosity function (LF) of the resolved sources are more consistent with 
sources found in other early- type galaxies than sources found in star- forming regions of galaxies. The 
source LF and the spectrum of the unresolved emission both indicate that there are a large number of 
unresolved point sources. We propose that these sources are associated with globular clusters (GCs), 
and that NGC 1600 has a large GC specific frequency. Observations of the GC population in NGC 1600 
would be very useful to test this prediction. Approximately 50-75% of the unresolved flux comes from 
diffuse gaseous emission. The spectral fits, hardness ratios, and X-ray surface brightness profile all 
point to two gas components. We interpret the soft inner component (a < 25", kT ^ 0.85 keV) as the 
interstellar medium of NGC 1600 and the hotter outer component (a > 25", kT ~ 1.5 keV) as the 
intragroup medium of the NGC 1600 group. The X-ray image shows several interesting structures. 
First, there is a central region of excess emission which is roughly cospatial with Ha and dust filaments 
immediately west of the center of NGC 1600. There appear to be holes in the X-ray emission to the 
north and south of the galaxy center which are roughly coincident with the lobes of the NGC 1600 
radio source. On larger scales, there is excess emission to the northeast, which we suggest may indicate 
the center of the group potential. The group galaxy NGC 1603 shows a tail of X-ray emission to its 
west which is probably due to ram-pressure stripping. 

Subject headings: galaxies: elliptical and lenticular — galaxies: ISM — intergalactic medium — 
X-rays: binaries — X-rays: galaxies — X-rays: ISM 



1. INTRODUCTION 

X-ray emission in early-type galaxies generally comes 
from two sources, hot {kT ^ 1 keV) interstellar gas, 
and hard X-ray point sources whose properties are 
consistent with low-mass X-ray binaries (LMXBs). 
X-ray bright galaxies (those with relatively high 
Lx/Lb ratios, where Lx is X-ray luminosity and 
Lb is blue optical lumino sity) are dominated by 
the interstellar gas (e.g., Forman, Jones. & Tucke 
119851 Trinchieri, Fabbiano, & Canizares 198 

while X-ray faint galaxies (low Lx/Lb) have 
a large proportion of emiss io n by LMXBs 
(iFabbiano Kim fc Trinchieril 119941: iPellegrinil Il994t 



llegr 

iKim et al."1996^ 'Sarazi n. Irwin, fc Bregmanli200(lD . 

bmce Fabbiano (198§7 

we have known that some of 

the off-nuclear X-ray point sources in spiral and ellip- 
tical galaxies have luminosities significantly exceeding 
the Eddington limit for 1 Mq . These ultra-luminous 
X-ray point sources (ULXs) appear to occur prefer- 
entially in star- forming regions, (e.g., the Antennae: 
iZezas fc Fabbiano! 2002); however, early-type galax- 
ies also contain bright point sources (e.g. . NGC 1399 
lAngelini. Loewenstein. fc Mushotzkvl l2001j) . In early- 



type galaxies, there appear to be X-ray point sources 
with Lx < 2 X 10'^^ erg s~^, consistent with accreting 
objec ts with masses < 15 Mq l|Irwin. Bregman. fc AthevI 
120041) . This mass limit is in line with current esti- 
mates of the upper mass limit of s tellar mass black 
holes for progenitor masses < 40 Mq ijFrver fc Kaloger'al 
2001). Above 2 x lO'^^ erg s~^, the number of sources in 
previously observed early-type galaxies may be consis- 
tent with the num ber of expected background sources 
ijirwin et al.l l2004|) . For this reason, we will adopt 
Lx = 2 X 10^^ erg s~^ as our minimum luminosity for 
a ULX candidate. Although the fainter ULXs could be 
explained by steady, spherically-symmetric, Eddington- 
limited accretion onto stellar-mass black holes, some 
other mechanism is required for the brighter ULXs. 

ASCA results have indicated that the total luminosity 
of LMXBs in early-type galaxies correlates better with 
the number of globular cl usters (GCs) than the opti- 
cal lu minosity of the galaxy ijWhite. Sarazin. fc Kulkarnil 
l2002^ . Chandra observations of early-type galaxies have 
also shown that a significant fraction (20%-70%) of the 
LMXBs are associated with globular clu sters i n the 
host galaxies fSarazin. Irwin, fc Bregman 2000, 20011 
Angelini et al.l i20Q_l; Kundu et al. 2003.' .Sarazin et aP 
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of the LMXBs were formed in GC s lIGrindlavl ITosl 
iSarazin et alJl20M IWhTte et alJl200l . and thus LMXBs 
can be used as tracers for GCs. 

In this paper, we discuss Chandra observations of 
NGC 1600, an X-ray bright E3 galaxy. NGC 1600 is the 
brightest member of the NGC 1600 group; NGC 1601 
{HQ away) and NGC 1603 (2f5 away) are the two near- 
est galaxies, b oth of which are non-interacting members 
(pc Vaucoulcu rs et al.ll992() . We adopted the distance t o 
NGC 1600 of 59.98 Mpc from lPrugniel fc SimienI |T996j) . 
which assumes H o = 75 km s~^ Mpc~-^ and uses the 
IFaber fc BursteiJ ljl988|) model that accounts for the 
Virgocentric flow and the Great Attractor. NGC 1600 is 
a boxy elliptical with a radially anisotropic, axisymmet- 
ric three-integral distribution function. Combined with 
its lack of significant rotation, its dynamics argue for a 
merger o rigin in which the effects of g as were not very im- 
portant f Matt hias fc Gerhardll999j) . iTerlevich fc Forbes! 
[2002) estimated the ag e of NGC 1600 to b e 7.3 Gyr 
with [Fe/H]= 0.41, while iTrager et alJ ^000) estimated 
an age of ~ 8.8 Gyr with [Fe/H]~ 0.15 through the inner 
5'.'7 and ~ 4.6 Gyr and [Fe/H]^ .24 through the in ner 
22'.'7. The colors (ISandagelll97a IFrogel et alJIT97l of 
NGC 1600 are consistent with colors from NGC 3379, a 
prototypical elliptical galaxy, suggesting star-formation 
has not occurred recently. In addition to X-ray emit- 
ting ga s, NGC 1600 also has cooler gas as indicated 
bv Ha llTrinchieri fc di Serego Ahghierilll99H) and dust 
ijFerrari et alJll99'^ Post-AGB stars seem capable of 
producing the necessary ionization/heating for the Ha 
and dust. 

In § 121 we discuss the observations and data reduc- 
tion of NGC 1600. After presenting the X-ray images 
in § 13 we discuss the properties of resolved sources in 
§ ^ The spatial distribution of the diffuse X-ray emis- 
sion and structures found in it are described in §|31 where 
they are compared to structures in other wavebands. We 
discuss the X-ray spectral properties of the sources and 
unresolved emission in § We estimate the gas and 
gravitational mass in § [T] Finally, we summarize our 
conclusions in §|S1 

2. OBSERVATION AND DATA REDUCTION 

NGC 1600 was observed in two intervals (observations 
4283 and 4371) on 2002 September 18-19 and Septem- 
ber 20 with live exposures of 26,783 and 26,752 s, re- 
spectively. The ACIS-35678 chips were operated at a 
temperature of —120 C with a frame time of 3.2 s. We 
determined the pointings so that the entire galaxy was 
located on the S3 chip, with the galaxy center offset from 
the node boundaries of the chip. Although a number of 
serendipitous sources were seen on the other chips, the 
analysis in this paper is based on data from the S3 chip 
alone. The data were telemetered and cleaned in Very- 
Faint mode, and only events with ASCA grades of 0, 2, 
3, 4, and 6 were included. Photon energies were deter- 
mined using the gain file acisD2000-08-12gainN0003.fits 
and corrected for time dependence of the gain^. We ex- 
cluded bad pixels, bad columns, and columns adjacent 
to bad columns or chip node boundaries. 

Chandra is known to encounter periods of high back- 
ground ( "background flares" ) , which especially affect the 

^ Seejhttp : //hea-www ■ harvard ■ edu/$\sim$alexey/ acis/ tgaln/ 1 



backside-illuminated SI and S3 chips^. We determined 
the background count rate from the SI chip to avoid the 
enhanced flux due to the galaxy on the S3 chip. Using 
Maxim Markevitch's lc_CLEAN progranP, we found the 
exposure intervals that were unaffected by background 
flares. The first observation showed clear evidence of a 
major flare in the first 20% of the observation. The sec- 
ond observation had some small fluctuations greater than 
20% from the mean rate. After these were filtered, ob- 
servations 4283 and 4371 had flare-free exposure times of 
21,562 and 23,616 s, respectively. We created a merged 
events file for imaging analysis after checking that the ob- 
servations were well registered; the separate events files 
were used for spectroscopic analysis. We took the back- 
grounds for extended regions from the deep blank sky 
backgrounds compiled by Maxim MarkevitcfP, adjusted 
them to the aspect histories of our observations, and 
changed the normalizations slightly to match the hard 
count rate (pha=2500:3000) of the blank-sky background 
with a relatively emission-free region on the S3 chip. For 
imaging analysis, we also included the "background" due 
to the readout artifact in ACIS using a script based on 
the MAKE_READOUT_BG progranP. 

We performed the data reduction and some of the data 
analysis using the Chandra analysis package CIAO 2.3'^ 
and NASA's ftools'*. Spectra were fit using xspec^, 
while correcting for the ACIS quantum efficiency (QE) 
degradation^ with the xspec ACISABS model. 

3. X-RAY IMAGE 

NGC 1600 has a combination of diffuse emission and 
resolved point sources. In order to image the diffuse emis- 
sion without deemphasizing the point sources, we adap- 
tively smoothed the background-subtracted, exposure- 
corrected image using a minimum signal-to-noise ratio 
(S/N) per smoothing beam of 3. Figure displays a 
true-color image of the adaptively smoothed X-ray im- 
age. This image was created by smoothing three ex- 
posure and background-corrected images in soft (0.3- 
1 keV), medium (1-2 keV), and and hard (2-6 keV) 
bands using the same kernel required for the total band 
adaptively smoothed image, and then combining them 
with the color coding red = soft, green = medium, and 
blue = hard. A logarithmic intensity scale was cho- 
sen to range between total band surface brightnesses of 
5x 10^^ counts arcsec"^ s^^ and 1 x 10^^ counts arcsec"^ 
s-i. 

In Figure 13 we show the corresponding Second Palo- 
mar Sky Survey (POSS II) optical (red) image using a 
linear gray scale. The field of view is the same as that 
in Figure ^ The largest galaxy near the upper cen- 
ter is NGC 1600, and corresponds to the brightest peak 
in the diffuse X-ray emission. Another group member, 
NGC 1603, lies to the east of NGC 1600. The galaxy 
north of NGC 1600 is the group member NGC 1601. On 
FigureEl the overlaid regions indicate the positions of the 
X-ray sources discussed in § 0] and listed in Table ^ All 
of the sources with fluxes determined to > 3cr (shown as 

See http : //cxc .harvard. edu/contrlb/ maxim/acisbg/| 
^ See http : //asc .harvard. edu/ciao2 . 3/ 

See http://heasarc.gsfc.nasa.gov/ftools/ 
^ See http : //heasarc . gsfc .nasa.gov/docs/sof tware/lheasoft/ 
^ See http://cxc.harvard.edu/cal/Acis/Cal_prods/qeDeg/ 
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Fig. 1. — Adaptively smoothed Chandra true-color S3 image 
(with red = 0.3-1 keV, green = 1-2 keV, and blue = 2-6 keV) 
of NGC 1600, corrected for exposure and background. The to- 
tal intensity scale is logarithmic and ranges from 5 x 10~^ counts 
arcsec"^ s~^ to 1 x 10""^ counts arcsec"^ s~^. The white square 
is the field of view of the Chandra S3 image. 
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Fig. 2. — Linear gray scale POSS-II Rod optical image of 
NGC 1600. The squares and circles indicate the positions of the 
X-ray sources with S/N> 3 and S/N< 3 in Table iTl respectively. 
Sources mentioned in the text are labeled. The dark square is the 
field of view of the Chandra S3 image. 



squares) are clearly seen in the adaptively smoothed X- 
ray image (Figure^, except for the central three sources 
that are embedded in strong diffuse emission. Very few 
of the weaker sources can be seen in Figure Q] 

The diffuse X-ray emission in Figure ^ shows several 
interesting structures. The central X-ray emission of 
NGC 1600 is elongated in a direction which is roughly 
aligned with the optical emission. However, there are 
some structures in the central X-ray emission, which are 
discussed below fS I5.1() . The emission around the ellipti- 
cal galaxy NGC 1603 in the east appears somewhat ex- 
tended, and there is a bridge of X-ray emission extending 
from NGC 1603 to the west toward NGC 1600. There 
is an X-ray source associated with the lenticular galaxy 
NGC 1601, but it is unclear whether it is extended in 
this image. On larger scales, there is a slightly elongated 
region of very extended X-ray surface brightness, with 
excess diffuse emission to the northeast of NGC 1600. 
(Adaptively smoothed images with the sources replaced 
by Poisson noise do not affect the gross morphology of 
the diffuse X-ray emission; that is, these diffuse features 
are not due to the smearing of point sources.) 

4. RESOLVED SOURCES 
4.1. Detections 
We used the wavelet detection algorithm (ciAO 
WAVDETECT program) with -\/2 scales ranging from 1 to 
32 pixels to identify the discrete X-ray source popula- 
tion on the ACIS S3. Since the wavelet source detection 
threshold was set at 10~^, <1 false source (due to a sta- 
tistical fluctuation in the background) is expected in the 
entire S3 image. Source detection was first performed on 
the separate observations to check their astrometric reg- 
istration; no significant offset was found. To maximize 
S/N, we analyzed the wavelet detection results from the 
combination of the two observations. We required our 
sources have wavdetect fluxes determined at the > 3(t 
level for all analyses except the identification of possible 



optical counterparts. 

At that level, the minimum detected count rate in 
the 0.3-6 keV band was 2.7 x 10"^ counts s~^; how- 
ever, the bright diffuse gaseous emission at the center 
of the galaxy makes it difficult to establish a minimum 
detectable flux over the entire image. For backgrounds 
of < 0.05 counts pix~^ and off-axis distances appro- 
priate to the S3 chip, sources with 20 counts should 
be detecte d at a roughly u niform completeness level 
(> 85%) l)Kim fc Fabbiand '20031 For sources satis- 
fying these criteria, the minimum detected flux was 
4.7 X 10~^ counts s'^. To determine source character- 
istics other than their flux, we used a local background 
with an area three times that of each source's wavdetect 
region. In a few cases of nearby sources, the source or 
background regions initially overlapped. We slightly al- 
tered these overlapping regions, preserving the ratio of 
areas and the net count rates. 

We also attempted detections in multiple bands (0.3-1, 
1-2, 2-6 keV) and compared detection rates to the total 
band (0.3-6 keV). In these sub-bands, 45, 68, and 34% of 
the total band sources were detected, respectively. There 
were two, one, and five sources detected in an individual 
sub-band that were not detected in the total band. None 
of these extra detections had fluxes which were significant 
at the > 3cr level. Two, eight, and five sources were 
detected in only one sub-band in addition to the total 
band, with two, six, and one of those sources having total 
band fluxes determined at > 3ct. In this X-ray bright 
galaxy, performing detections by sub-bands provided no 
advantage. 

In Table ^ we list all discrete sources detected by 
WAVDETECT ovcr the 0.3-6 keV range. The sources are 
ordered by increasing projected radial distance d from 
the center of the galaxy. Columns 1-8 provide the source 
number, the lAU name, the source position (J2000), the 
projected radial distance and semi-major distance a from 
the center NGC 1600, the wavdetect count rate with its 
Icr error, and the S/N for the count rate. The fluxes 
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were corrected for exposure and the instrument PSF. The 
first three sources are clearly extended; when the count 
rate was determined using I'.'S circular regions centered 
on the source centroid, none of those sources are signif- 
icant at the 3ct level. Although the position of Source 
1 is close to the optical/IR nucleus of NGC 1600, we 
are not confident it is a point source. Therefore, we 
adopted the 2MASS Point Source Catalog position of 
R.A. = 4'^31'"39!87 and Dec. = -5°5'10'.'5 as the location 
of the center of NGC 1600. All listed positions include 
astrometry corrections based on optical/IR counterparts 
(S I4.2|I : the overall absolute astrometric errors are proba- 
bly ~0'.'5 near the field center, with larger errors further 
out. 

In addition to the <1 false source in the entire S3 
field of view, some of the detected sources may be 
foreground or (more likely) background objects unre- 
lated to NGC 1600. If we consider the fluxes of all 
our detected sources, we expe ct ~48 unrela t ed so urces 
based on the source counts in ' Brandt et alJ l)2000D and 
iMu s hotzky e t al. (2000). Using the minimum detected 
> 3(7 flux, we expect ~15 unrelated sources. Incom- 
pleteness will reduce both of these estimates, especially 
the first number. Unrelated sources should be spread 
out fairly uniformly over the S3 image (Figure^, except 
for the reduced sensitivity at the center due to bright 
diffuse emission and at the outer edges of the field due 
to reduced exposure and increased PSF. Sources close 
to NGC 1600 are more likely to be associated with the 
galaxy, while sources far from NGC 1600 are more likely 
unrelated to NGC 1600. 

4.2. Identifications 

Sources in Table H were cross-correlated against op- 
tical/IR catalogs to identify possible counterparts and 
to improve the absolute astrometry of the observa - 
tions. We used the Tycho-2 Catalog ijHog et alJl200(]D . 
2MASS'' Point Source and Ext ended Source C atalogs^, 
and the USNO-B Catalog (,Monet et al.ll2003j) to iden- 
tify seven optical/IR counterparts to the X-ray sources. 
Four of these sources were used to determine the as- 
trometry: Source 18 corresponds to USNO-Bl 0849- 
0044132, an i? = 19.4 mag object with a non-stellar 
PSF; Source 23 corresponds to NGC 1601, a nearby 
lenticular galaxy north of NGC 1600; Source 33 corre- 
sponds to 2MASS 04313613-0503081, a J = 16.8 mag 
star; and Source 50 corresponds to Tycho-2 4742-254-1, 
a B = 11.7 mag star. After correcting for an astromet- 
ric shift of about 0'.'5, the residual astrometric errors are 
ft! 0'.'5. 

We did not use the remaining three sources for astrom- 
etry since they were associated with extended X-ray or 
optical emission. The central source (Source 1) is » 0'.'9 
from the optical center of the galaxy (after the astro- 
metric correction discussed above). As the X-ray source 
is extended, we used a 1'.'5 radius circular region cen- 
tered on the X-ray source to measure or limit the fiux 
of any central point source; this gave a count rate of 
(7.6 ± 2.9) X 10^^ counts s~^ (Icr error bars). Since this 
source was neither particularly hard nor detected very 



significantly as a point source, we conservatively adopt 
its 3(7 upper limit luminosity of 6.7 x 10^^ erg s~^ as an 
upper limit to a central active galactic nucleus (AGN). 

Source 40 is « 1'.'7 from the 2MASS center of 
NGC 1603, a nearby elhptical galaxy east of NGC 1600. 
There is clearly extended X-ray emission centered east 
of the source position that appears more coincident with 
the galaxy center. Since NGC 1603 extended well be- 
yond this X-ray source (20 mag arcsec"^ isophotal K 
fiducial elliptical aperture semi-major axis, r_k20fe, of 
18'.'4), Source 40 may be a point source offset from the 
center of NGC 1603. Similarly, Source 63 is « 2'.'9 from 
2MASX 04312667-0507309, a galaxy candidate whose 
r_k20fe is 7'.' 2; Source 63 may also be a point source offset 
from this galaxy's center. Although Source 10 is close to 
an optical point source in Figure |21 the X-ray source is 
« 3'.'4 south of the optical object; thus, we do not con- 
sider this to be a reliable optical identification for the 
X-ray source. 

4.3. X-ray Luminosities and Luminosity Functions 

To convert the source count rates into unabsorbed X- 
ray luminosities, we used the adopted Chandra X-ray 
spectrum 16.11 Table [21 row 3) of the resolved sources 
to convert 0.3-6 keV count rates into 0.3-10 keV flux. 
We then assumed each source was at the distance of 
NGC 1600, 59.98 Mpc, yielding a conversion factor of 
4.12 X 10^2 erg count-^ Column 8 of Table[T]hsts the 
X-ray luminosities in units of 10^^ erg s"""^, which range 
roughly from 2.8 x 10^^ to 4.8 x 10*° erg s"^ Since 
Sources 33 and 50 are likely foreground stars, their lumi- 
nosities are probably overestimates. 

By restricting the sample to the sources with a uni- 
form completeness of 85% (> 20 net counts implying a 
count rate limit > 4.7 x 10~* counts s~^, and d > 40" 
corresponding to a background < 0.05 counts pix~^, 
ii l4.1|l . excluding the sources corresponding to NGC 1601 
(Source 23) and NGC 1603 (Source 40), and excluding 
the very bright source corresponding to a foreground Ty- 
cho star (Source 50), we created our analysis sample of 
20 sources (see Notes in Table Pi . We expect 11 ±2 fore- 
gr ound/background so urce s based on the source co unts 
in [Brandt et alJ lj2000jl and lMushotzkv et alJ l|200(]|) . 

In Figurc|21 we display the cumulative luminosity func- 
tion (LF) of our analysis sample. The LF should be the 
sum of the point source (LMXB) population of NGC 1600 
and the foreground/background population. We fit the 
LF using the same techniq u es we have used previously 
llSarazin et all l2?M l2001t |l rwin. Sarazin. Sz BregmanI 
120021) : a single power-law, broken power-law, and a cutoff 
power-law were all used to model the LMXB population. 
The background source population was modeled as dis- 
cussed in the previous references. A single power-law fits 
the data very well based on the Kolmogoroff-Smirnov 
(KS) test: 

dN 



NoL 



See jhttp : //www ■ Ipac . caltec h.edu/2inass/releases/second/doc 

* When a source appeared in both 2MASS catalogs, the Point 
Source Catalog positions were used. 



where L39 = Lx/(10'^^ erg s~^). The best fit was 
determined by the maximum likelihood method, and 
the errors (90% confidence interval) were determined by 
Mnntp Carlo techniques. We found iVo = 21.llJ'n Q and 

/explsup.html^^ . , , , , . „ i • ' j- 

uL — ^.uu_Q 35. Although a cutoff power-law is not re- 
quired by the fits, the total luminosity of all of the sources 
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Fig. 3. — Histogram of the observed cumulative luminosity 
function of resolved sources in our analysis sample. The contin- 
uous curve is the sum of the best-fit LMXB luminosity functions 
(eq. and the expected background source counts. The vertical 
line indicates the completeness limit of our sample. 



diverges at the high luminosity hmit for the best-fit LF. 
The lack of very bright sources (beyond those observed) 
in NGC 1600 can be explained by Poisson fluctuations; 
however, it seems likely that the correct underlying LF 
of the sources in NGC 1600 is either a bit steeper than 
the best-fit value, or that it steepens or has a cut-off at 
high luminosities beyond those observed. 

Early-type galaxies tend to have broken or cut- 
off power-law LFs with the break or cutoff oc- 
curring well below the luminosities measured in 
NGC 1600, while star forming galaxies tend to 
have single power-law LFs that extend to low lu- 
minosities ISarazin et al. 2001; Blanton et al. 2001; 
Finoguenov & Jones 2002; Zezas & Fabbiano 2002; 
Jeltema et al.. .2003; S ivakoff. Sarazin. fc Irwini i2003i: 
Randall. Sarazin. fc Irw in 2004). The high luminosities 
of the sources in NGC 1600 make it difficult to directly 
compare its luminosity function to other galaxies. 
Therefore, one must either extrapolate the NGC 1600 
LF down to lower luminosities or extrapolate the 
LF of other galaxies to higher luminosities in order 
to compare them. The differential LF slopes at the 
highest luminosities in early-type galaxies tend to be 
steeper than 2.5, and tend to cluster around 1.5 in 
galaxies with some level of star formation. The best-fit 
slope of NGC 1600 is intermediate between the two; 
however, within the errors, the slope is more consistent 
with the other early-type galaxies. The normalization 
of NGC 1600 is higher by at least a factor of four 
than in any of the galaxies in the above references. 
Extrapolating the best-fit luminosity function, we would 
find ~ 40 and ~ 420 sources above 5 x 10^^ erg s~^ 
and 5 x 10^^ erg s~^, typical break luminosities and 
minimum observed luminosities in other early-type 
galaxy observations. If there is a break in the luminosity 
function below 2 x 10'^^ erg s^^, the numbers of fainter 
sources would be reduced. 

Only one of the sources in the analysis sample is fainter 
than 2 x 10'^^ erg s^^. There are six additional sources 
which could qualify as ULX candidates. Sources 1-3, 7, 
52, and 71. Since Sources 1-3 appear extended, we used a 



I'.'S radius circular region to measure the flux of a possible 
point source. The count rates of Sources 1 and 2 are suffi- 
cient to qualify as ULX candidates; however, their fluxes 
are not determined at the 3cr limit. Of the remaining 
candidates, only Sources 7 and 71 are ULX candidates 
with well-determined fluxes. Source 7 is not included in 
the analysis sample since it is only 35'.'2 from the nucleus 
of NGC 1600. Source 71, a source at the edge of the chip 
did not make it into the analysis sample due to its small 
number of counts, ~ 17; its luminosity is boosted by a 
large exposure correction. Including Sources 7 and 71, 
we observe 21 ULX candidates. At the flux limit corre- 
sponding to Lx = 2 x 10'^^ erg s~^, we expect 11 ±2 fore- 
gr ound/backgrqund so urces based on the source co unts 
in iBrandt etall 1)2000(1 and iMushotzkv et all l|2000 ^ for 
the entire chip. The number of ULX candidates in ex- 
cess to the expected background is 10 ± 5. Even if we 
consider sources with Lx > 4 x 10'^^ erg s~^, there 
are ten ULX candidates and only w5 unrelated fore- 
ground/background sources are expected. This corre- 
sponds to an excess of 5 ± 3.3^ The error budgets of 
both excesses are dominated by Poisson counting errors 
in NGC 1600. 

In a sample of 28 early-type galaxies, the num- 
ber of sources with Lx > 2 x 10^^ erg s~^ (0.3- 
10 keV ) was equivalent to th e expected numbe r of fore- 
ground/background sources ijirwin "eTaLllmil i". Al- 
though the number of ULX candidates is greater than 
the average number of unrelated sources, the number 
of ULX candidates in NGC 1600 could be a result of 
cosmic variance in the foreground/background sources. 
One item in support of this possibility is that the fit- 
ted slope of the luminosity function is consistent with 
the typical sl opes of foreg round /background lu minos- 
ity functions l|Brandt et alJl200H iGiacconi et alJl2001(l . 
To examine the possibility of cosmic variance, we com- 
pared the source densities we observed to the 0.5-2.0 keV 
source densities of the Chandra Deep Fields""^^. 

Within two D25 (a < 147'.'3), the observed source den- 
sity of ULX candidates in NGC 1600 is 2250±796 deg-^, 
while the remaining field of the S3 chip has a source 
density of 884 ± 245 deg^^. In order to compare the 
background source densities with t hose in IBrandt et alJ 
(|2001) and IGiacconi et all pOOl . we converted our 
count rates to energy fluxes in the 0.5-2.0 keV band 
using a r = 1.4 power-law spectrum. Then, the mini- 
mum detectable 0.5-2.0 keV flux for our observation of 
NGC 1600 is 1.36 x 10" erg cm'^ s"^ The Chandra 
Deep Fields predict foreground/background source den- 
sities of - 690 ± 260 deg-2 (North) and 520 ± 110 deg-^ 
(South). Assuming a background source density of 
600 ± 100 deg~^, there is an ~ 2.1tT excess within two 
D25 of the galaxy corresponding to an excess of ~ 6 ± 3 
sources. The errors are dominated by counting errors in 
NGC 1600 as opposed to uncertainty in the background 
source density. The excess near the galaxy favors an as- 

^ We use \/N statistics since they better represent the true low- 
count Poisson lower limit Icr confidence level than the Gehrels ap- 
proximation of 1 -I- V-'V + 0.75 (Compare with eq. [7], eq. [11], and 
Table 2 in|Gehrcls 11986)). 

This study used surface brightness profil e distances which a re 
consistent with Ho = 74 ± 4 km s"! Mpc'^ ITonrv et aTll2001^ . 

We chose to use the 0.5-2.0 keV band since all but one of the 
ULX candidates were detected in soft as well as hard bands. 
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sociation of ULX candidates with NGC 1600; however, 
cosmic variance cannot be ruled out. 

When one subdivides the ULX candidates into a fainter 
sample Lx = 2-4 x 10^^ erg s~^ and a brighter sample, 
Lx > 4 X 10'^^ erg s~^, we find that the brighter sam- 
ple source densities are more consistent with the back- 
ground than the fainter sample. In particular, there 
are no sources in the bright sample, which has better 
completeness than the faint sample, within two effec- 
tive radii; in the faint sample, there are four sources. 
One possible explanation is that the overabundance of 
the bright ULX candidates is due to cosmic variance 
of foreground/background sources, while the overabun- 
dance of the faint ULX candidates is from sources within 
NGC1600. The ULX candidates could then be brought 
more in line with the findings of Irwin et al. (;2p04) by a 
coincidence of cosmic variance and a 40% overestimate 
in the distance to NGC 1600. Such large differences do 
occasionally occur between recessional velocity distances 
and surface brightness fluctuation distances; however, in 
the case of NGC 1600, this would require NGC 1600 have 
a peculiar velocity of ~ +13GGkm s^^. Since this pecu- 
liar velocity is more typical of cluster infall and there is 
no nearby large cluster, we do not think this model is 
likely. 

If cosmic variance does not explain the excess number 
of bright sources, then these bright sources are actually 
ULX candidates associated with NGC 1600. NGC 1399 
has three sources with Lx > 2 x 10'^^ erg s~^, the 
largest number at these luminosities among previously 
observed early-type galaxies. NGC 1407 may have five 
such sources: how ever, its di stance is highly unce rtain 
prwin et al.ll2003(l . Although Ueltema et all l|2003D find 
that NGC 720 has nine ULX candidates, only one has a 
luminosity > 2 x 10"^^ erg s~^ for Hq = 75 km s~^ Mpc~^. 
Thus, we believe that NGC 1600 may have the largest 
number of ULX candidates brighter than 2 x 10^^ erg s^^ 
observed in an early- type galaxy to date. 

There are a variety of models explaining the pres- 
ence of ULX sources. Some of these include favor- 
able viewing angles of anisotropic radiation (King et al. 
l200l . super-Eddington accretion by high mass X-ray bi- 
naries (HMXBs) at the thermal-timescale mass transfer 
rate feing 20 03), accretion onto in termediate-mass black 
holes (iColbert fc Mushotzkvlll999|) . super-Eddington ac- 
cretion of LMXBs in the soft X-ray transient state 
IIKingll2002D . and super-Eddington accretion fro m thin 
accre tion disks of stellar mass black holes (,Begelmanl 
12002(1 . If anisotropic radiation was the cause of ULXs 
in NGC 1600, one would predict the existence of a large 
population of ULXs seen at fainter fluxes, due to mis- 
alignment with the preferred axis of the anisotropic radi- 
ation, in addition to the intrinsically faint and isotropic 
radiating LMXB populations; our observation doesn't go 
deep enough to observe such sources. 

Since NGC 1600 has an observed age of 7.3 Gyr 
( Terlevich & Forbes 2002) a nd its photometric colors 
I Sandage 1973; Frogel et al] 1197^ are typical for an 
elliptical galaxy, one would not expect its X-ray bi- 
naries to be HMXBs. Comparisons of the FIR de- 
tections of NGC 1600 from IR AS in the 60m and 
100m bands with the radio flux IfBirkinshaw fc DaviesI 
((1985') 1 could be consistent with either AGN or stella r 
heating of interstellar dust ijCondon fc Brodericklll99H) . 



iTrinchieri fc di Serego Alighieril l|1991[l and lFerrari et alJ 
l|1999l) both found that post-AGB stars are the likely 
source of the ionization of the gas and heating of the dust 
in NGC 1600. Thus, the FIR detections do not necessar- 
ily indicate that there is recent star formation as needed 
for HMXBs. Additionally, the cumulative spectrum of 
the ULX candidates in NGC 1600 is not consistent with 
the disk blackbody model found to fit w ell in ULXs as- 
sociated with spiral galaxies ( Makishima et al.ll2000|) . 

Intermediate- mass black holes may be created if the 
proge nitor mass of a star is > 40 Mq l|Frver fc Kalogerl] 
120011) and that star sinks to the center of a globu- 
lar cluster where it can grow up to ~ IO^Mq in 10^" 
years (Miller fc Hamilton 2002). Under this model, 
NGC 1600 would be expected to have a large number of 
globular cluste rs. The LMXB-GC conn ection in early- 
type galaxies llSarazin et alJ 120001 1200 It j Angelin i et alJ 
l200lUKundu et al.ll2003l: iSarazin et alJl2003i|) means that 
LMXBs in the soft X-ray transient state or black hole 
LMXBs with thin accretion disks would also be consis- 
tent with a large number of globular clusters. Thus, it 
seems plausible that the large population of ULXs in 
NGC 1600, if they are not due to cosmic variance in the 
number of unrelated sources, requires a large population 
of GCs. Unfortunately, the globular cluster population 
of NGC 1600 does not appear to have been determined. 

4.4. Hardness Ratios 

Hardness ratios or X-ray colors are useful for crudely 
characterizing the spectral properties of sources, and can 
be applied to sources that are too faint for detailed spec- 
tral analysis. We determined the observed X-ray hard- 
ness ratios for the sources, using the same tec hniques we 
used previously ( Sarazin et al. 2000, 2001: Bla nton et alJ 
rwin et alJI2002() . We define three hardness ratios 
as H21 = {M-S)/iM + S), H31 = {H-S)/{H + S), and 
H32 = {H - M)/{H + M), where S*, M, and H are the 
total counts in the soft (0.3-1 keV), medium (1-2 keV), 
and hard (2-6 keV) bands, respectively. As compared to 
our previous definitions, we have reduced the hard band 
from 2-10 keV to 2-6 keV. Since the 6-10 keV range is 
dominated by background photons for most sources, this 
should increase the S/N of the hardness ratio techniques. 
The hardness ratios measure observed counts, which are 
affected by Galactic absorption and QE degradation in 
the Chandra ACIS detectors. In order to compare to 
other galaxies, it is useful to correct the hardness ratios 
for these two soft X-ray absorption effects. Therefore, we 
have calculated the intrinsic hardness ratios, denoted by 
a superscript 0, using a correction factor in each band ap- 
propriate to the best-fit spectrum of the resolved sources. 
The intrinsic hardness ratios and their Ict errors are listed 
in columns 10-12 of TableHl 

Although we have plotted H31 vs. H21 in the past, 
plots of H32 vs. H21 allow for better separation of sim- 
ple spectral models. In Figure 01 we plot both H31'' vs. 
H21'' and H32° vs. H21° for the 20 sources in the anal- 
ysis sample. The hardness ratios for the sum of those 
sources are (H21°,H320,H3lO) = (-0.26,-0.39,-0.59); 
the uncorrected hardness ratios are (H21,H32,H31)= 
(-hO.il, -0.33, -0.24). Sources with ~40 net counts had 
errors similar to the median of the uncertainties, ~ 0.2. 
The errors scale roughly with the inverse square root of 
the net counts. 
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Fig. 4. — Hardness ratios for the sources in our analysis sample. 
Here, H2l0 = (M° - S°)/{M° + S°), mi° = {H° - S°)/{H° + 3°), 
and H320 = (H" - M0)/{H° + M^), where S° , M", and are 
the counts in the soft (0.3-1 keV), medium (1-2 keV), and hard 
(2-6 keV) bands, corrected for the effect of Galactic absorption 
and QE degradation according to the best-fit spectra of resolved 
sources. The area of each circle is proportional to the observed 
number of net counts. The solid curve and large diamonds show 
the hardness ratios for power-law spectral models; the dashed curve 
and small diamonds show the ratios for intrinsic absorption of 
Njj = 4 X 10'^^ cm"'^; the diamonds indicate values of the power- 
law photon number index of F = (upper right) to 3.2 (lower left) 
in increments of 0.4. Both models underwent the same correction 
as the sources. The Ic error bars at the upper left illustrate the 
median of the uncertainties. 



In previously studied galaxies, most of the sources lie 
along a broad diagonal swath extending roughly from 
(H21,H31) - (-0.3,-0.7) to (0.4,0.5). Usually, these 
hardness ratios were not corrected for Galactic absorp- 
tion and QE degradation; the latter effect was not known 
at the time of some of the previous studies. Since there 
are fewer sources and the absorption/degradation cor- 
rections tend to push sources to the lower-left part of 
the diagram, this swath is less evident in NGC 1600. In 
Figure 0] the solid line corresponds to hardness ratios 
for power-law source spectra with F =0-3.2. In calcu- 
lating these model hardness ratios, Galactic absorption 
and QE degradation were applied to the model spec- 
tra, and the hardness ratios were corrected for these 
effects using the best-fit spectral model of the resolved 
sources as described above. The dashed line corresponds 
to a similar model, with an intrinsic absorbing column 
of 4 X 10^^ cm~^. Most of the sources are consistent 
with F =1.6-2.2, with the sum of these sources corre- 
sponding to F « 1.7. The majority of the sources lie 
roughly between the two models in H31° vs. H21°. In 
1132*^ vs. H21°, we see that most of the sources have an 

H32'' 0.5, but the two models still roughly contain 

the sources. When compared to the 1131" vs. 1121" plot, 
the effects of absorption in the II32" vs. H2l" plot are 
larger, and the effects of absorption and varying power- 
law index are more nearly orthogonal. Although it is 
suggestive that the sources tend to lie on one of the two 
tracks, the errors are large and the number of the sources 
is small, so this could be a coincidence. 

In NGC 4697 and the bulge of NGC 1291, a total of 
four sources had (II21,H31) « (—1,-1). Scaling from 
the number of sources in NGC 4697, one would have only 
expected ~2 supersoft sources in NGC 1600, so the lack 
of any strong candidates is not surprising. Moreover, the 
supersoft sources in NGC 4697 would all have been below 
our detection limit at the distance of NGC 1600. Another 
problem is that the soft X-ray response of the Chandra 
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ACIS S3 detector was much worse at the time NGC 1600 
was observed due to QE degradatioiP. Although the 
hardness ratios are corrected for this effect, the correction 
is based on the average (hard) spectrum of the sources; 
this correction would be too small for supersoft sources. 
Supersoft sources in NGC 1600 would have been hard to 
detect unless they were very bright. 

Among the sources with > 3tT known fluxes, we find 
three sources (42, 58, and 60) with very hard spec- 
tra, (H21",H31") > (0.5,0.5), and three sources (17, 
36, and 50) with little hard emission, H31° ~ —1.0. 
The very hard sources may be unrelated, strongly ab- 
sorbed AGNs, similar to the sources which produce the 
hard component of the X-ray background, and which 
appear strongly at the faint fluxes in t he deep Chan- 
dra observations of bl ank fields (.Brandt et al.l l2000t 
Mushotzk v et all 120001: iGiacconi et all 120011) : however, 
all three sources have large hardness ratio errors. The 
sources without hard emission may also be unrelated 
foreground/background sources as indicated by studies 
of other galaxies ( „ Sara zin et al. 2001) and deep blank 
s ky i mages (e.g., IGiacconi et al.l [20011) : Source 50 (See 
S I4.2|I is clearly a foreground object, while Source 36 has 
large hardness ratio errors. 

The hardness ratios of Sources 1-3 are softer than a 
majority of the sources in the analysis sample; however 
the errors of Sources 2 and 3 could bring them more 
in line with the analysis sample. Since these extended 
sources appear softer, it is more likely that they are small 
scale structures in the diffuse gas (See § I5.2|l , than that 
they are a number of confused sources. 

4.5. Variability 

With two observations we could test for variability on 
two time scales, intra-observation, and inter-observation. 
Since the mid-points of the observations were only 60 
ks apart, the timescales are not remarkably different. 
We searched for variability in the X-ray emission of 
the resolved sources over the duration of the Chandra 
observations individually and jointly using the KS test 
(see Sarazin ct al. 2001). Additionally, we compared the 
count rates between the two observations. In Tabled we 
report the sources that had a > 95% variability proba- 
bility in any one of these four tests. Among the 43 3cr 
flux detected sources away from the edge of the chip. 
Sources 1, 8, 24, 30, and 59 were apparently variable. 
Source 1 had a 2.1cr (96% probability) rate decline in 
count rate between the two observations; however, the 
variability was not significantly detected within the in- 
dividual observations. This suggests that there may be 
a point source variable on a timescale > 20 ks buried in 
the diffuse emission of this extended source. The flux in 
a 1'.'5 circular region centered on the position of Source 
1, which includes the galaxy center, is nearly constant 
between the two observations. This suggests that the 
buried variable point source, if it exists, is not a central 
AGN. For Source 8, the KS test on the joint observa- 
tions yielded a 95% probability of variability. Source 24 
shows a variability at the 99% probability in the first 
observation; in this observation, six of its seven counts 
occur within ~ 4 ks. Source 30 undergoes a 2. Oct (95% 
probability) rate increase between the two observations. 
Finally, Source 59 has a 98% probability of being variable 
in the joint observation KS test. With 43 sources, ~ 2 
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false positives are expected at a > 95% limit. Of the four 
sources variable at this limit, Sources 24 and 59 show the 
strongest behavior of variability, both in the form of out- 
bursts. Although Source 71 shows 95% variability in the 
first observation, its lightcurve is irregularly sampled due 
to its location at the edge of the chip. Source 13, a weakly 
detected source, had a very high variability probability, 
99%, in the second observation since all three counts in 
that observation occurred in the last 750 seconds. The 
remaining weak sources that are marked as potentially 
variable did not exhibit such clearly variable behavior as 
in Source 13. 

5. UNRESOLVED X-RAY EMISSION 

As mentioned in § El the morphology of the unresolved 
X-ray emission is complex. This emission is a combina- 
tion of unresolved point sources and diffuse gaseous emis- 
sion. Ideally, we would like to use the spectral properties 
of the gas and resolved sourc e s (See [i El to disentan- 



^le th e two as in lSarazin et al.l (|2001') and'Sivako ff et aP 
()2003tl : however, the changing diffuse gas temperature 
of NGC 1600 does not allow this separation. For this 
reason, we chose to use the entire 0.3-6 keV band in 
analyzing the unresolved emission of NGC 1600. From 
the spectral fits, we estimate unresolved point sources 
contribute 30% of the total unresolved counts. For 
comparisons with previous studies, we left in resolved 
sources with S/N< 3. As Source 1 may be the peak of 
the unresolved emission, we did not exclude it. Since 
Sources 2 and 3 were extended, we only excluded regions 
corresponding to 1'.'5 circular regions located at their cen- 
ters. We used the blank-sky background to statistically 
remove background events. 

5.1. Radial Profile of the Unresolved X-ray Emission 

In order to compare the spatial distribution of the un- 
resolved X-ray emission in NGC 1600 with the optical 
emission by stars, we adopted the Third Reference Cata- 
logue of Bright Galaxies (RC3) values for the optical pho- 
tometry's effective radius {ves = 45'.'4), position angle 
{PA — 15), and ellipti city (e = 0.324), wh ich assumes a 
de Vaucouleurs profile ijde Vaucouleurs et al. 1992). The 
corresponding s emi-m ajor axis, OoS is 55'.'2. Although 
iRembold et al.l l)2002f ) find a different profile in the in- 
frared (a Sersic profile with n ~ 1-2, an effective radius 
of ~10-15", a semi-major axis position angle of ~10-20, 
and ellipticity of 0.17-0.34), the values of the position 
angle and ellipticity are roughly consistent with the op- 
tical profile. The corresponding semi-major axis is 17'.'6 
using the optical ellipticity. Since the effective radius is 
only used for scaling in this paper, the large discrepancy 
between the two radii does not make a difference in our 
study. 

We determined the surface brightness profile (SEP; 
Figure EJ and the hardness ratios (Figure for a se- 
ries of elliptical annuli with semi-major widths of 5", ex- 
tending to 180". We could have used much narrower an- 
nuli near the center of NGC 1600, but there is significant 
non-axially-symmetric structure in the image there fS I5.2l 
below). In the SEP, the dotted line displays the (opti- 
cal) best-fit de Vaucouleurs profile with effective semi- 
major axis fixed at 55'.'2 and the dashed line displays the 
(J-band) best-fit Sersic n = 1.65 profile with effective 
semi-major axis fixed at 17'.'6. The normalizations of the 



optical profiles were varied to achieve the best fit. It is 
clear that neither the optical nor the infrared profile fit 
the unresolved X-ray emission well. 

We first tried to fit the SEP using a single beta model 
profile, 

Ix{a) = Io ' (2) 

where Oc is the core semi-major axis. However, that 
model was clearly rejected with a of 431 for 33 degrees 
of freedom (dof). The best-fit single beta model has a 
small core radius a^, leading to essentially a power-law 
model. The model underestimates the emission in the 
inner 20" and beyond 90" while overestimating the emis- 
sion between 20" and 60". Figure |S1 clearly shows that 
there are at least two components to the X-ray SEP of 
NGC 1600. Thus, we tried to fit the SEP with a double 
beta model: 

2-\ -3A„„er-(-i/2 

a 

1 



Ix{a) 



^0, inner 



1 



^c, inner 



Q-c, outer 



-3/3outcr + l/2 



(3) 



The core radius of the outer component Oc, outer is poorly 
determined and consistent with zero. That is, the 
outer part of the SEP could be fit by a power-law sur- 
face brightness; however, a pure power-law form would 
produce a peak in the X-ray SEP at the very center 
which is not seen in the image. Therefore, we chose 
to freeze the outer core semi-major axis at 25". We 
then found /o, inner = 0.97tQ;Qy counts arcmin^^ 



a 



c, inner — -^^■'^ — 2.3 ' 



arcsec (4.2 kpc in projection), /Jjimer = 
l.mto.li ^.outer = 0.039t^:™^ counts s-1 arcmin-2, 
Oc, outer = 25.0 arcscc (7.3 kpc in projection), and 
Pouter = 0-36to oi- The values of this fit were consistent 
with the error bars from the double beta model profile 
with a free outer core semi-major axis, and the of this 
fit, 39.0 for 30 dof, was only 0.9 higher. Therefore, we 
adopted this fit, the solid-line in Figure El as the best-fit 
SEP model. To ensure that the inclusion of Source 1 did 
not affect the fit, we attempted to fit the same model 
without using the first annulus. The new fit was within 
the error bars of the above fit. 

A (3 oi 1.18 is much larger t h an in most X-ray bright 
galaxies ijForman et alJ 119851 iTrinchieri et alJ Il986'l . 
However, these previous fits were done with a single beta 
model and on much poorer resolution data. Eelow, we 
will argue that the inner beta model component of the 
SEP is interstellar gas in NGC 1600, while the outer 
beta model component is intragroup gas. It may be that 
the intragroup medium has compressed the galactic gas, 
decreasing the value of Oc and steepening the profile (in- 
creasing /?). 

Figure El shows the profiles of the hardness ratios for 
the unresolved emission. As was true for the SEP, the 
hardness ratio profiles show a break in behavior at a ~ 
25". The inner regions of the galaxy are noticeably softer 
than outer regions. It is unlikely that this trend is due 
to unresolved point stellar sources. The optical profile 
of the galaxy is more centrally condensed than the X- 
ray SEP (Figure O, so unresolved stellar sources should 
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Fig. 5. — Surface brightness profiles, with la error bars, of the 
unresolved emission (0.3-6 keV) as a function of projected semi- 
major axis a. The dotted and dashed curves show the RC3 de 
Vaucouleurs profile and J-band Sersic profile, respectively, with 
fixed effective radii and normalizations varied to fit the X-ray sur- 
face brightness. The solid curve is the best-fit double beta model 
profile. All fits were for a < 180". 



contribute more in the inner regions. Yet, the average 
spectrum of the resolved sources is much harder than any 
of the diffuse emission, and would not cause the central 
regions of the diffuse emission to be softer. A positive 
temperature gradient and/or abundance gradient in the 
diffuse gas in the inner regions is the likely cause for the 
softer emission; in S 16.21 below, we show that the X-ray 
spectra do indicate that this change in the hardness is 
due to a temperature gradient. 

5.2. Structural Features in the Unresolved X-ray 
Emission 

In Figure [7| we display images where the resolved 
sources were replaced by the appropriate local Pois- 
son noise. In each image, a cross marks the center 
of NGC 1600. On the left, we display the adaptively 
smoothed image of the total X-ray emission (0.3-6 keV) 
with logarithmic gray scale ranging from 1 x 10~^ counts 
arcsec"^ s~^ to 2.9x 10"** counts arcsec"^ s~^. From this 
image, we subtracted the double beta elliptical SBP; the 
excess emission image is shown in the middle panel with 
both limits of the gray scale reduced by a factor of two. 
The tail extending from NGC 1603 to the west toward 
NGC 1600 is shown clearly. This image shows that there 
is extended, diffuse emission around both of the smaller 
galaxies NGC 1601 and NGC 1603. There is some inter- 
esting residual structure near the center of NGC 1600, 
which is discussed in more detail below (§ 15. 3() . The 
structure near the center includes an excess emission re- 
gion to the west of the center of NGC 1600. On larger 
scales, there is excess X-ray emission to the east and 
northeast of NGC 1600. In the right panel of Figure[71 we 
display the hardness ratio image (5" Gaussian smoothed) 
of H41 = {H + M-S)/{H + M + S) between -0.4 (black) 
and 0.4 (white). This displays the softer central emission 
well. One can also see that the NGC 1603 tail appears 
softer than its surroundings. These features are also ev- 
ident in Figure n 

To quantify the surface brightness differences associ- 
ated with the features in Figure we determined the 
X-ray surface brightness of the unresolved emission in a 
number of elliptical annular pie regions. For each annu- 
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Fig. 6. — Hardness ratios with Itr error bars for the unresolved 
emission as a function of semi-major axis a. Here, H21 = (M — 
S)/{M+S), H31 = {H-S)/{H+S), and H32 = (H-M)/(H+M), 
where S, M, and H are the observed counts in the soft (0.3-1 keV), 
medium (1-2 keV), and hard (2-6 keV) bands. 



lus, we subdivided the region by PA (Figure|Hl). To quan- 
tify non-uniform morphology, we first fit a constant line 
to the surface brightnesses at each annulus (the dashed- 
line), iteratively excluding PA regions from the fit that 
diverged by more than 2a. In the innermost annulus, 
emission from unresolved point sources and any inaccu- 
racies from using the optical isophotes will be largest. 
Additionally, there appear to be small depressions in the 
emission just to the north and the south of the galaxy 
center, and a larger depression at the east edge of the 
annulus. These make it difficult to accurately establish 
a baseline for this annulus. There is a clear excess of 
emission [PA = 255° - 345°, > 4.3ct) to the west of the 
second annulus. This excess continues in the third an- 
nulus {PA = 255° — 315°, > 3. Oct), where an excess also 
appears 180 deg away {PA — 105 — 135, a > 5.6). In the 
fourth and fifth annulus, the western excess has disap- 
peared, while the eastern excess begins to take up most 
of the northeast quadrant {PA = 30° - 120°, > 4.0ct). In 
the fourth annulus, there is also an isolated excess in the 
south (~ 3.5ct). Finally, the last annulus clearly shows 
the NGC 1603 tail (- 4.5cr). 

The excess to the northeast of NGC 1600 in annuli 
3-5 of Figure |H1 may be due to gas associated with the 
NGC 1600 group potential rather than the galaxy. We 
note that this northeast excess starts at about the same 
radius where the X-ray surface brightness profile in Fig- 
ure has an inflection point and the larger scale com- 
ponent dominates. This suggests the gas at these radii 
is responding to a potential with a larger scale. Both of 
these suggest that the outer gas is bound to the group, 
and that the center of the group potential is to the north- 
east of NGC 1600. 

We believe that the tail west of NGC 1603 is the re- 
sult of ram pressure stripping. NGC 1603 is at a pro- 
jected distance from NGC 1600 of ~44.5 kpc with its 
velocity redshifted from NGC 1600 by - 284 km s'^. 
Since NGC 1603 is part of the NGC 1600 group, we 
assumed that they were at the same distance. Given 
the tail seen in the X-ray image, it is unlikely that the 
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Fig. 7. — (Left) Adaptively smoothed Chandra S3 image (0.3-6 keV) of NGC 1600, with sources removed and corrected for exposure 
and background. The gray scale is logarithmic and ranges from 1 X 10~^ counts arcsec"'^ s~^ to 2.9 X 10~* counts arcsec"^ s~^. (Middle) 
The image on the left minus the best-fit double beta model surface brightness profile. The gray scale is also logarithmic and ranges from 
5 X 10~^ counts arcsec~2 s~i to 1.4 X 10""* counts arcsec"^ s~i. (Right) The image of the hardness ratio //41 = (H + M — S)/{H + M + S), 
Gaussian smoothed with ct = 5". The gray scale is linear and ranges from -0.4 (black) to 0.4 (white). In each image, a cross marks the 
center of NGC 1600. 



only component of the velocity of NGC 1603 relative to 
NGC 1600 is along the line of sight. We assume that 
the two transverse components of the relative velocity 
are each about the same as the line of sight component. 
Thus, the total relative velocity of NGC 1603 is about 
« a/3 X 284 km s"^ « 490 km s"^ This would be a rea- 
sonable value for a circular orbital velocity for NGC 1603 
around NGC 1600, whose observed ra dial velocity disper- 
sion IS 321 km s-i (|Faber et al.l lT98^. We argued above 
that there may be a significant potential associated with 
the NGC 1600 group, which would increase the estimated 
velocity. 

We estimated the gas densities and pressures in 
NGC 1603 and in its environment to see if ram pressure 
could be sufficient to strip gas and form the tail. For 
the gas in NGC 1603, we assumed three uniform density 
spherical annuli with widths of 2" each. Since the X- 
ray emission from NGC 1603 is softer than the center of 
NGC 1600, we assume a spectrum with kT = 0.6 keV and 
solar abundances for the gas in this galaxy. With these 
assumptions, we find that the electron number densities 
of the gas in NGC 1603 are 4.9 x 10-^ 1.9 x lO'^, and 
5.4 X 10~^ cm~^ and the pressures, Pgas, are 9.4 x 10~^\ 
3.6 X 10"", and 1.0 x 10"" dyne cm^^ in the 0-2, 2-4, 
and 4-6" annuli. We estimated the density in the group 
gas around NGC 1603 from the X-ray surface brightness 
in a hemispherical annulus centered on NGC 1603 (radius 
from NGC 1600 = 6"-12", PA = 0°-180°). This surface 
brightness agrees within the errors with the prediction by 
the modeled surface brightness profile at the semi-major 
distance of NGC 1603 from NGC 1600. We assumed 
a spectrum with kT = 1.5 keV and solar abundances. 
This gave an ambient gas electron density of 4.0 x 10~^ 
cm~^ at the projected semi-major distance of NGC 1603. 
Assuming a relative velocity of 490 km s^^, the ram pres- 
sure of the ambient gas would be Pram ~ 1.6x10^^^ dyne 
cm~^. The ram pressure would be lower if NGC 1603 
is at a larger radius than its projected radius (i.e., in 
front of or behind the center of the NGC 1600 group). 
Roughly speaking, the condition for ram pressure to strip 
the gas is that Pram > Pgas, which appears plausible in 
the outer 4-6" of NGC 1603. Since the tail points to 



the west towards the higher density gas near NGC 1600, 
the ambient gas density may have been higher in the 
past, making ram-pressure stripping easier than it is at 
its present position. 

The NGC 1603 tail has a softer spectrum (Figure [7| 
right) than the surrounding gas. Hardness ratios from 
counts in the same pie annuli used to construct Figure |H1 
confirm that the tail is softer in H21 by « 2.6ct. Thus, the 
jump in emission due to the tail is dominated by a jump 
in soft emission, consistent with ram-pressure stripping 
of cooler galactic gas by hotter group gas. 

NGC 1601 also has a velocity redshifted from 
NGC 1600 by ~ 309 km s^^, and its projected distance 
from NGC 1600 is only ^28.5 kpc. Since the emission 
associated with NGC 1601 is fainter, its estimated ther- 
mal pressure is smaller than in NGC 1603. This, com- 
bined with its slightly higher velocity, suggests that ram- 
pressure stripping should be stronger in NGC 1601 com- 
pared to NGC 1603. Although there is a hint of a tail 
toward the east of NGC 1601 in Figure [3 it has only a 
1.3ct significance when the surface brightness in a hemi- 
spherical annulus east of NGC 1601 is compared to a 
similar annulus west of NGC 1601. There are a number 
of possible explanations of why we do not see a significant 
tail around NGC 1601 while we see one for NGC 1603. 
First, NGC 1601 might have a small transverse veloc- 
ity; the tail would then be projected onto NGC 1601 
along our line of sight. Additionally, the tail could be 
too faint, the projected velocity could be smaller than 
we estimated for NGC 1603, or projection effects could 
mean NGC 1601 is actually in a lower density region in 
the group gas than we calculate from the surface bright- 
ness profile and its projected radius. 

5.3. Central X-ray Structure and Multi-Wavelength 
Comparisons 

The X-ray emission near the center of the galaxy is 
elongated similarly to the optical emission (Figuresdand 
121). However, the details of the central X-ray structure 
are complex. Figures [7| and |S1 display the excess emission 
immediately west of the galaxy center, as well as holes in 
the emission immediately north and south of the galaxy 
center. 
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One of the purposes of this Chandra observation was to 
compare the X-ray structure to the extended emission- 
line filaments and dust. In Figure IHl we display the 
Ha -|-[N II] image (Trinchicri & di Scrcgo Alighicri 199^) 
overlaid on the excess emission image of the inner 140 x 
140". The peak of the Ha corresponds to excess emission 
west of the galaxy center; however, the detailed struc- 
ture of the Ha and excess X-ray emission do not have a 
one-to-one correspondence. IMacchetto et all |l996) find 
'--^ 5 X 5" of Ha-|-[N II] emission centered on the optical 
center, with a slight north-south elongation. There is 
excess X-ray emission centered on the galaxy; however, 
a correlation of excess emission is made difficult by the 
possible presence of a central AGN and holes in the X- 
ray emission ^ 2'.'5 to the north and south of the galaxy 
center. 

On the left of Figure^l we display the excess emission 
image of the inner 4 x 40". It is overlaid b y 4.885 GHz 
radio contours from iBirkinshaw fc Davie3 |^85). The 
double radio lobe structure is oriented in the north-south 
direction, and the lobe positions are roughly consistent 
with holes in the X-ray emission. Although the south- 
ern lobe's centroid is a bit offset from the hole, the 5" 
Gaussian smoothed image of H41 indicates that the lobe 
is cospatial with a region of softer emission. Considering 
the smoothing scale, this could be consistent with cool 
gas surrounding a radio bubble. No variation in H41 is 
observed near the northern lobe. These holes in X-ray 
emission may be due to the radio lobes, although this is 
uncertain. 

In Figure 1101 we also display a.V — R reddening map 
(iFerrari et al The darker areas of the color in- 

dex map correspond to regions of larger V — R and ex- 
tinction Av- There is a clear filamentary structure to 
the west corresponding to a mean Ay = 0.034 ± 0.030 
(jFcr rari et al. 1999), where the standard deviation is 
due to non-uniform extinction over the area of mea- 
surement. The optical extinction can be converted 
into a hydrogen c olumn densit y assuming Nh = 5.9 x 
W^Av/Rv cm-^ llSnitzeHlTpTl . Taking i^y = Av/{B- 
V) = 3.2, there is an excess Nr of 0.63±0.55 x 10^° cm"^ 
beyond the Galactic value of 4.86 x 10^*^ cm~^. If one 
assumes that the X-ray emission of a mekal model with 
T = 0.85 keV and solar abundance is absorbed by this 
extra column, the surface brightness would be reduced 
by only compared to that expected with no excess 
absorption. Of course, we actually observe an excess in 
the X-ray emission in this region of slightly larger ex- 
tinction. The excess emission to the west of the galaxy 
is ~ 1.5-2. 5 times the values expected from the surface 
brightness profile model. 

Both the dust lanes and emission line filament appear 
to be cospatial in projection with the enhanced X-ray 
emission. There are several possible explanations for 
such a correlation. It might be that the hot and cool 
gas are in thermal contact, and that hea t is being con - 
ducted from the hot gas to the cool gas fSDarkg''199^. 
This might cool the X-ray gas; if it remained at nearly 
the same pressure due to the weight and pressure of the 
surrounding hot gas, the density and X-ray emissivity of 
the X-ray gas would increase. Mixing between the hot 
and cool gas could have a similar effect, as long as the 
X-ray gas didn't mix to a cool temperature out of the 
X-ray band. Third, the cool gas might result from radia- 
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Fig. 8. — Surface brightness profiles, with Icr error bars, of 
the unresolved emission (0.3-6 keV) as a function of PA for an- 
nuli with varying projected semi-major radii a. The dashed line 
indicates the best-fit constant surface brightnesses of each annulus, 
iteratively excluding PA regions more than 2a from it. The last 
annuli excludes PA regions that extend beyond the S3 chip. 




Fig . 9. — Ho -|-[N II] contours ITrinchieri fc di Serego Alighieril 
11991]) overlaid on adaptively smoothed image of unresolved X-ray 
emission (140 X 140") with a best-fit elliptical model for the surface 
brightness profile removed. The contours indicate line emission of 
0, 1.7, 2, 4.1, 6.2, 8.3, 12.4, 18.6, 27.6, 41.4 xlO"!'' erg cm'^ 
arcsec"'^. The cross indicates the center of the galaxy. 




Fig. 10. — Left, adaptively smoothed gray scale image of un- 
resolved X-ray emission (40 X 40") with the best-fit beta-|-power- 
law surface brightness profile removed. Contours are radio (4.885 
GH z: 1, 2, 3, 4, 5, 6 mjy) with the beam in the lower-right cor- 
ner IBirkinshaw fc Davieijlll985h . The cross indicates the center of 
the galaxy. Right, gray scale image o{ V — R (0.625, 0.63, 0.635, 
0.64 ); the darker (redd er in color) values in correspond to larger 
Av IFerrari et aTlraggh . 
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tive cooling of X-ray gas. The cooling timescale at the 
center of NGC 1600 is on the order of 300 Myr, which 
is consistent with this explanation. On the other hand, 
it might be difhcult to understand the presence of dust 
grains in cooled X-ray gas. 

6. SPECTRAL ANALYSIS 

We extracted spectra of the resolved sources and 
diffuse emission in NGC 1600, restricting analysis to 
the 0.7-9 keV range. The lower limit was taken to 
avoid calibration uncertainties, while there are few non- 
background counts beyond 9 keV. Since the telescope 
collecting area changes very rapidly near 2 keV and cali- 
bration problems led to poor fits in that spectral region, 
we chose to excise the 1.9-2.1 keV band when we fit line 
models (mekal). We do not find that this edge signifi- 
cantly affected fits of continuum models. All of the spec- 
tra were grouped to have at least 25 counts per spectral 
bin prior to background correction to enable our use of 
statistics. In some cases, the fitted range of the spec- 
tra did not extend up to 9 keV because there were too 
few counts to form a bin up to this limiting energy. 

The results of the spectral fits are summarized in Ta- 
ble [3 Spectra were extracted for the resolved point 
sources ('Sources') and the unresolved diffuse emission 
excluding the point sources ('Unresolved'). The third 
column gives the geometric region for the spectrum; 
'Field' implies the entire S3 chip. The value of the 
absorbing column density (Nh) applied to all compo- 
nents of the model emission spectrum is given in col- 
umn 4. In this and other columns, values in parenthe- 
ses are fixed (not allowed to vary ). The fixed value 
of Nh is the Galactic value from IDickev fc Lockm'ai] 
((1990). To correct for the QE degradation, we used 
the XSPEC ACISABS model. Other elliptical galaxies are 
known to have both h ard (point sou rce) and soft (dif- 
fuse gas) components ijSarazin et al.ll2()0l|) . Under the 
hard component, we fit three different models, thermal 
bremsstrahlung ('bremss'), power-law ('power'), and disk 
blackbody ('diskbb'). We always used the mekal model 
for the emission spectrum from hot diffuse gas. For the 
hard component, columns 5-7 give the spectral model(s), 
the temperature (for bremsstrahlung or disk black- 
body) or photon number spectral index F, and the unab- 
sorbed fiux of the hard component(s), Fi^ (0.3-10 keV). 
Similarly, columns 8-10 give the temperature T^, overall 
heavy element abundance relative to solar, and flux for 
the soft mekal component (s). For the unresolved emis- 
sion, the spectra exclude regions around each of the > 3a 
flux determined resolved sources. The last two columns 
give the number of net counts in each spectrum, and 
per dof for the best-fit model. All errors reported in the 
spectral analysis are 90% confidence level errors. Brack- 
ets are used when either the upper or lower bound on the 
confidence interval was unconstrained. 

The background spectrum for the resolved sources were 
determined locally, using the same nearby regions as dis- 
cussed in § 0] For the spectra of the unresolved emis- 
sion, we used the deen blank sky backgrounds compiled 
by Maxim MarkevitclP. 

The spectra of several spatial regions were analyzed. 
We have adopted the ellipticity and PA of the optical de 
Vaucouleurs fit. Some of the regions are scaled to the el- 
liptical optical isophote containing one-half of the optical 



light, the "one effective radius" region. In NGC 1600, the 
semi-major axis of this isophote, Ocfi, is 5 5'.' 2. Since there 
were few source counts near the center of the galaxy, 
we fit the sources for the entire field. We searched for 
changes in unresolved emission with radius using annuli 
containing ^ 1000 net counts. 

6.1. X-ray Spectrum of Resolved Sources 

All resolved source spectral analysis was performed 
on sources with > 3a determined fluxes. Sources 1- 
3 were excluded because we are unsure that they are 
point sources. Sources 23 (NGC 1601), 40 (NGC 1603), 
and Source 50 (CSC 04742-00254) were also excluded 
because they are clearly not associated with NGC 1600. 
Since there are few source counts near the galaxy's cen- 
ter, we first extracted the cumulative spectrum of re- 
solved sources within three effective radii, yielding only 
372 net counts. This was insufficient to produce a well- 
constrained fit, so we chose to fit the sources in the en- 
tire field. The observed spectrum, containing 1318 net 
counts, is shown in Figure ITTI 

The combined spectrum of the sources was reasonably 
well-fit by either a thermal bremsstrahlung model with 
kTfi = 4.73 keV (Table 12 row 1) or power-law model 
with a photon number spectral index of F — 1.76 (row 
3). The fits were not improved significantly when the 
absorbing column was allowed to vary (rows 2 and 4), 
so we fixed the hydrogen column at the Galactic value, 
Nh = 4.86X lO^o cm^^ ( Dickcv & Lockman 1990). Since 
the sources were very luminous, we also attempted a disk 
blackbody model; however, its was worse by ~ 20 for 
the same number of dof. Finally, we attempted a disk 
blackbody -I- power-law spectral model {kTdisk = 1.44, 
r = 1.95, power-law responsible for ~ 70% of flux). This 
model had a lower than our bremsstrahlung or power- 
law fits by ^ 2; however, the dof was also reduced by 
2. Additionally, this fit had an unconstrained disk tem- 
perature, and poorly constrained power-law. We could 
have adopted either the bremsstrahlung or the power-law 
model for our best-fit. For comparison to other papers 
and since the power-law was constrained more tightly 
than the bremsstrahlung temperature, we adopted the 
power-law model as our best-fit. This model, includ- 
ing Galactic absorption and QE degradation, is shown in 
FigurelTTlwith the residuals to the fit. The power-law in- 
dex is consistent with what was found from the hardness 
ratios. It is softer than the best-fit value of the sources 
simultaneously fit in a survey of 15 early- type galaxies 
(|Irwin et al.ll200:^ : however, it is consistent with both 
the spread of indices found when the galaxies are fit sep- 
arately and the softer indices of the higher luminosity 
subsamples. 

Approximately half of the sources in the entire field 
are expected to be foreground/background sources. On 
the other hand only '--^ 1 is expected to be a fore- 
ground/background source within three effective radii. 
To ensure that unrelated sources do not heavily bias the 
spectral fits, we compared fits to the less constrained 
three effective radii spectrum with the fits to the entire 
field spectrum. Since all of the fits were consistent within 
the errors of the entire field spectrum fits, we do not be- 
lieve the foreground/background sources heavily bias the 
spectral fits. 
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Fig. 11. — Upper panel: the cumulative X-ray spectrum of 
the resolved sources with > Sit detected flux in the entire field of 
NGC 1600. The spectrum has la error bars and is overlaid by the 
solid histogram of the best-fit model spectrum (Table l2l. Lower 
panel: the contribution to with the sign indicating the sign of 
the residual. 



6.2. X-ray Spectra of Unresolved Emission 
6.2.1. Projected Spectra 

For the unresolved emission, the spectrum of the inner 
effective radius had 3662 net counts. This spectrum is 
shown in FigureEl First, we attempted to model the un- 
resolved emission with a soft mekal component represent- 
ing the emission by diffuse interstellar gas (Tabled row 
11). The was large, 168.6 for 87 dof. Since the unre- 
solved emission includes unresolved point sources as well 
as diffuse gas, we added the adopted best-fit spectrum 
of the resolved sources to model the unresolved sources. 
The fit was much improved, = 110.0 for 86 dof; how- 
ever, it was still rejected at the > 95% level. We allowed 
the absorbing column to vary, but this did not signifi- 
cantly improve the fit and was still consistent with the 
Galactic value. Therefore, we have assumed Galactic ab- 
sorption for the remaining fits. Since the hardness ratios 
of the unresolved emission indicated spectral evolution 
at a ~ 25", we attempted a two-temperature gas solu- 
tion (row 9). We found a good fit, = 78.7 for 83 dof; 
however, the abundance of the low temperature gas was 
unconstrained. Therefore, we tried a two-temperature 
gas model with the abundances tied together. This fit, 
with a low mekal temperature of 0.85 ± 0.04 keV, a high 
mekal temperature of 2.55^Q gg keV, and an abundance of 
1.071 

a40 solar, was almost as good as when both abun- 
dances were free. We also attempted a mekal cooling 
flow model; however, its fit, = 94.4 for 85 dof, was 
worse than the two-temperature model. Therefore, we 
adopted the two-temperature gas model (row 10) as our 
best-fit model for unresolved emission (Figure I12|l . Al- 
though the best-fit flux of this model from the unresolved 
point sources was not large, the upper limit on the flux 
of 1.41 X 10~^^ erg cm~^ s~^ in 0.3-10 keV band stiU in- 
dicates unresolved sources could be a non-trivial source 
of emission in the inner effective radius. 

Since this galaxy was bright enough to flt multiple an- 
nuli of unresolved emission, and there were indications of 
multiple temperatures of the diffuse gas, we attempted 
to determine the radial dependence of gas temperature 
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Fig. 12. — Upper panel, the cumulative X-ray spectrum of the 
unresolved emission within 1 a^g of NGC 1600. The notation is 
the same as in Figure [TTI 



and abundance. We used eight annuli each with approx- 
imately 1000 net counts, out to a = 180". In Table [2 
rows 11-18, we show the results of assuming the model 
was the sum of the resolved point source model with its 
normalization free and a single temperature gas mekal 
model. Most of the fits produced reasonable x^- Any ra- 
dial changes in the abundance are dwarfed by the errors 
of the fits; although the best-fit abundances were mainly 
subsolar. On the other hand, it is clear that the first 
two annuli have a much lower temperature, 0.85 keV, 
than the outer five annuli, ^ 1.5 keV. The temperature of 
the third annulus is intermediate between the two tem- 
peratures. These fits are in rough agreement with the 
two-temperature model found within one effective radius, 
a < 55'.'2; however, the flux of the unresolved sources 
is smaller and the temperature fit for the hotter gas is 
larger when fitting the inner effective radius at once, as 
opposed to in multiple annuli. The best-fit normaliza- 
tion of the power-law and mekal model suggest that for 
a < 180" they both contribute approximately equally to 
the flux. Within looff, the diffuse gas is dominant by at 
least two-to-one. From the best-fit fluxes of the annular 
flts, we can estimate the X-ray luminosity (0.3-10 keV) 
to be ^ 2.9 X 10^^ erg s~^ in gas and ^ 2.4 x 10"*^ erg s~^ 
in unresolved sources for a < 180". Around 30-40% of 
the gaseous luminosity comes from the cooler gas. The 
resolved 3cr sources in the entire field have an X-ray lu- 
minosity of ~ 1.4 X lO''^ erg s~^. In a < 180", this scales 
to ^ 0.4 X 10^^ erg s^^. Bolometric corrections increase 
the source luminosity by 44% and the total gaseous lu- 
minosity by 28%. 

Using the RC3 optical profile, a < 180" corresponds 
to ~ 80% of the tota l optical light, L B,m% = 8.5 x 10^° 
Lbq l|0'Sullivan. Forbes, fc Ponmar]l200l|) . The source 
X-ray-to-optical ratio is 3.3 x 10'^° erg s~^ L^p-,! approxi- 
mately four times th at found in NGC 46 97 ( Sarazi n~et alJ 
120011) or NGC 1553 ('Bl anton et alJ200ll|) . It is more than 
ten times t he expected contr i bution from discrete sources 
reported in lO'SuUivan et ani)200lD . Since the unresolved 
sources dominate the total source luminosity, there is an 
excess in the flux of faint sources in addition to the excess 
in the number of detected sources. It is unlikely that cos- 
mic variance would lead to an excess number of sources 
at all fluxes. This suggests that the number of LMXBs 
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found in NGC 1600 is not proportional to the number of 
stars as estimated by the optical light, or that the un- 
resolved source flux is overestimated. The most likely 
solution is that NGC 1600 has a high specific frequency 
of globular clusters; this is discussed in more detail later 
in the paper (See §|H1). 

The large unresolved hard component flux in the spec- 
tral models suggests that there are many unresolved 
point sources. This is also roughly consistent with the 
observed luminosity function for resolved sources fii l4.3|) . 
If we extend the best-fit luminosity function down to 
10^^ erg and adjust the normalizations to account 
for the observed sources within Sooff, we expect an un- 
resolved fiux of ~ 2.2 X 10"^'^ erg cm^^ s~^ in the 0.3- 
10 keV band. Although this is approximately half of 
the fiux estimated by the spectral fits, it is within the 
errors of the spectral fit flux. This agreement between 
the luminosity function and spectral fits could point to a 
large number of unresolved point sources. It might also 
suggest that there is no strong break in the luminosity 
function at lower X-ray luminosities than observed. 

6.2.2. Deprojected Spectra 

We also fit deprojected spectra. We did the depro- 
jections by fitting the spectra from the outside annulus 
to the inside. Each annulus was fit assuming the best 
fit model for the emission in each of the outer annuli, 
and assuming simple geometric projection to include the 
emission from outer annuli in the inner ones. First, we 
allowed the power-law normalization and mekal model 
to vary within each annulus. The outer four annuli had 
unphysically large abundances, and in general the abun- 
dances were poorly constrained. The temperature results 
were similar to those for the projected fits. Next, we fit 
the deprojected spectra assuming a solar abundance in 
all annuli to ensure that the abundance problems in the 
previous deprojection did not cause problems in the tem- 
perature determination. This deprojection also produced 
results similar to the projected fits, except the jump be- 
tween a low temperature gas and a high temperature gas 
moved inward. Since these results are qualitatively the 
same as the projected fit, but we had to assume a con- 
stant abundance, we choose not to present the details of 
these results. 

6.2.3. Galaxy and Group Gas 

Both the unresolved emission within one effective 
radius and the unresolved emission in annuli out to 
180" point to a model involving two components of 
hot gas. In the inner 25-40" region, there is gas 
with kT ^ 0.85 ke V. With a stellar velocity disper- 
sion of 321 km s-i l|Faber et alllTosl . NGC 1600 has 
a stellar kinetic temperature of ~ 0.65 keV, predict- 
ing an X-ray t e mpera t ure of ^ 1 zfc 0.2 keV ba sed on 
iDavis fc Whiig IpOl : TBrown fc Bresfmanl l(T99l . This 
prediction is consistent with the measured X-ray tem- 
perature. At semi-major distances beyond 40", the gas 
is hotter, kT ~ 1.5 keV. This temperature is more con- 
sistent with the gas tempe rature of X-ray bright groups 
(jHelsdon fc Po nman"2 000l) than X-r ay luminous early- 
type galaxies fO 'Sulliv an et al.ll2003D . Since the abun- 
dances are poorly constrained, we cannot tell if there is 
any gradient in the abundance, or if the average abun- 
dance is higher or lower than solar. 



Near the same radius where there is a transition in 
spectrum, we saw transitions in the behavior of the sur- 
face brightness profile and the hardness ratios. This 
break in behavior of the surface brightness, hardness ra- 
tio, and spectrum at a ~ 25-40" (7.3-11.7 kpc in projec- 
tion) might be explained if emission from a lower temper- 
ature diffuse gas from the interstellar medium (ISM) of 
NGC 1600 dominates in the inner regions and a higher 
temperature diffuse gas from the intergalactic medium 
(IGM) of the NGC 1600 group provides the emission at 
larger radii. 

The total gas X-ray-to-optical ratio is 3.4 x 
2^q30 gj.g g-i ■ the galaxy gas luminosity is ~30-40% 
of the total gas lu minosity. Compared to th e Lx,hoi-L b 
relationships from O'SuUiv an et al.| l)2001|) . the galaxy 
and the total gaseous luminosity are consistent with the 
standard deviation around the best-fit relationship to 
the early-type galaxies excluding AGNs, brightest cluster 
galaxies and dwarfs. The total emission is consistent with 
the best-fit relationship for the brightest group galaxies. 

7. MASS DETERMINATION 

By combining the SBP model and the spectral fits of 
the diffuse gas, we estimated the gas and gravitational 
mass around NGC 1600. In eq. JSJ, the SBP is repre- 
sented as the sum of two functions. Because the X-ray 
emissivity is a quadratic function of the density, eq. Q 
does not lead to a simple expression for the gas density 
as the sum of two terms. Thus, we assumed that the two 
terms in the SBP represented two cospatial phases of gas, 
a 0.85 keV gas following the inner beta model profile and 
a 1.5 keV gas following the outer beta model profile. To 
remove the ~ 30% contribution to the counts by unre- 
solved point sources (See § EJ, we multiplied the SBP 
by 0.7. We integrated the sum of the physical densities 
over volume to determine the gas mass, and we used hy- 
drostatic equilibrium to estimate the gravitational mass. 
For both estimates, we assumed the galaxy was shaped 
as an oblate spheroid. 

In Figure lT^ we display the resulting gravitational and 
gas mass profiles. At a ~ 20", there is a slope change in 
the gravitational mass profile. This indicates the likely 
presence of two gravitational potentials, supporting the 
view that the exterior gas is group gas and the interior 
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Fig. 13. — The estimated gas and gravitational mass profiles 
around NGC 1600. The dotted line and right axis display the gas 
mass profile, while the solid line and left axis display the gravita- 
tional mass profile. 
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gas is galaxy gas. On the other hand, the total gravita- 
tional mass within a < 180" is only 2.0 ± 0.5 x lO^^M© 
with errors in the hot gas temperature dominating the 
error budget. The mass-to-light ratio for a < 180" is 
M/Lb = (24 ±6) Mq/Lbq. Within one effective radius, 
the mass-to-light ratio is (10 ± 3) Mq/ Lbq- The total 
mass is not that much larger than might be expected for 
NGC 1600 alone, which may not be consistent with the 
suggestion that the outer component is due to a group 
dark matter potential. 

The total gas mass within a < 180" is 1.14 ± 
0.06 X 1O^°M0, with errors in the surface brightness 
dominating the error budget. Assuming a normal el- 
liptical galaxy stellar mass loss rate per stellar lumi- 
nosity of 1.3 x 10~"Mq/Lbq (jMathews fc Brighentl 
1200 3'). the rate of stellar mass loss for the entire 
galaxy is about 1.4 M© yr~^. Almost all of the gas 
mass around NGC 1600 can be attributed to stellar 
mass loss in NGC 1600, assuming an age of 7.3 Gyr 
l|Terlevich fc Fo rbes"2002}) . This is not surprising since 
NGC 1600 dominates the stellar population of the 
NGC 1600 group. On the other hand, it is not con- 
sistent with a large amount of intergalactic group gas in 
the system. 

8. CONCLUSIONS 

We have used Chandra observations to study the X- 
ray emission from the point sources and unresolved emis- 
sion of NGC 1600. Since NGC 1600 is - 60 Mpc away, 
we could only resolve the very brightest point sources. 
We detected 71 total point sources, of which 45 were 
bright enough to have fluxes determined at > 3a. We 
identified two of the sources with foreground stars, two 
with known galaxy members of the NGC 1600 group, 
and three with non-stellar objects, presumably galaxies. 
We did not clearly detect a point source at the center of 
NGC 1600, and we put a conservative 0.3-10 keV upper 
flux limit of 6.7 X 10'^^ erg s~^ on the luminosity of the 
central AGN. 

Twenty-one of the sources without clear associa- 
tions to objects besides NGC 1600 are brighter than 
2 X 10^9 erg s"! (0.3-10 keV); approximately 11 ± 2 
of those sources are expected to be unrelated fore- 
ground/background sources. The excess is strongest 
within two D25, where we observe a source density 
of 2250 ± 796 deg~^ compared to the expected 600 ± 
100 deg~^. NGC 1600 may have the largest number of 
ULX candidates observed to date in an early-type galaxy, 
although cosmic variance in the background source popu- 
lation cannot be completely ruled out. A combination of 
cosmic variance for the high flux sources and a distance 
overestimate of ~ 40% could also reduce the number of 
ULX candidates; however, this model requires NGC 1600 
have a peculiar ve locitv of ^ -|-1300km s~^. As found in 
llrwin et al.l l)2003j) . the bright LMXBs have softer spectra 
than is typical for fainter LMXBs in early- types galaxies. 
The spectrum of sources in NGC 1600, T = 1.761^5;^, 
and the luminosity function slope, a = 2.00^J'35, both 
suggest that these sources are different than ULX candi- 
dates observed in star-forming galaxies. The age of the 
galaxy and the X-ray spectra of the sources argue against 
these sources being HMXBs. 

Since this galaxy is X-ray bright, we expected the un- 
resolved emission to dominate the flux. Although this 



is true, the contribution from unresolved point sources 
is nearly as large as that due to the diffuse gas in the 
best-fit models. Even if one uses the lower limit on the 
unresolved point source flux, unresolved point sources 
are responsible for ^25% of the unresolved flux. Com- 
bined with the large normalization of the X-ray luminos- 
ity function, this suggests that NGC 1600 has a larger 
population of LMXBs tha n found in m ost prev iously 
obser ved galaxies (Sarazin et alJ l2001t Blanton et alJ 
2001; Finoeuenov & Jones 2002; Zczas & Fabbiano 2002 
Jeltem a ct al. 2003; Sivakoff et al. 2003; RandaU et all 
2003)- When normalized to the optical Imirinosity of 
the galaxy, the source X-ray-to-optical ratio is also very 
high, Lx/Lb = 3.3 x 10^° erg s"^ L~q. This ratio is a 
factor of two larger than found in most early-type galax- 
ies. Recently, Gilfanov (2004) suggested that total X-ray 
source luminosities correlate better with near-infrared lu- 
minosities than with optical luminosities. From the ex- 
trapolated 2MASS Ks magnitude (8.04), Mk,,q = 3.39, 
and the distance of NGC 1600, we find that the source 
X-ray-to-near-infrared ratio is almost ari_ order of magni- 
tude higher tha n the average found in|Gilfanov (2004). 
iGilfanovl l|200l limited the sources he considered in his 
total X-ray luminosities to Lx > 10^^ erg s~^, so part 
of the difference might be explained if NGC 1600 has a 
very large number of fainter sources. 

Previous studies have shown that a large fraction of 
LMXBs in elliptical galaxi es are located in glo bular 
dusters USarazin et a,].il200nl l2n0lHAnge]ini et aIj l200Tt 
iKundu et alJ l2003t iSarazin et alJ 120031) . and that the 
luminosity of LMXBs may correlate bett er with the 
GC population than the o ptical luminosity l) White et alJ 
l2002t ISarazin et aI]l2003D . Unfortunately, we were un- 
able to find a determination of the GC population in 
NGC 1600, which is more distant than most galaxies 
which have GC observations. However, one explanation 
for the large number of LMXBs in NGC 1600 would be 
that this galaxy has a large GC specific frequency. Thus, 
we predict that NGC 1600 will be found to have large 
numbers of GCs. Observations of the GC population in 
NGC 1600 would be very useful to test this prediction. 
Such a large population of GCs might be connected to 
NGC 1600's position as the central elliptical galaxy in a 
group with a significant potential; it may be more closely 
related to cD galaxies, which have larger GC populations 
(e.g., Harris 1991) than normal ellipticals. 

We note that many of the bright X-ray sources asso- 
ciated with NGC 1600 are at large distances from the 
galaxy center. One possibility is that these sources are 
due to cosmic variance in the background source popula- 
tion. For the sources within two D25, the variance would 
have to be much larger than the ~ 8% cos mic variance 
betw een the Chandra Deep Fields found by ijRosati et alJ 
l2002^ . Another possibility is that these sources are actu- 
ally located in intergalactic GCs in the NGC 1600 group, 
rather than being directly associated with NGC 1600. 
Recently, intergalactic GCs have probably been detected 
in nearby groups and clusters ('e.g.. .Bassino et al. 2003) . 

The X-ray image, surface brightness profile, and spa- 
tially resolved spectra suggest that there are two compo- 
nents to the gas around NGC 1600. There is a component 
which is centered on NGC 1600, which has a small spatial 
scale (< 25") and a lower temperature (0.85 keV), which 
we propose is gas which is bound to NGC 1600. A second 
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component is centered to the northeast of NGC 1600, has 
a larger spatial scale, and is hotter (1.5 keV); we argue 
that this gas is bound to the dark matter potential of 
the NGC 1600 group. The X-ray image suggests that 
the center of the potential of the group is slightly dis- 
placed from the center of NGC 1600. One would expect 
that NGC 1600 would be moving in the group poten- 
tial, and this motion could affect the distribution of the 
galaxy and group gas. 

We also observe structure to the diffuse emission on 
small scales. Excess emission is seen west of the center 
of NGC 1600 in the inner 35". This excess emission is 
partially cospatial with Ha and dust filaments. Possible 
models to explain this correlation of hot and cooler gas 
include thermal conduction between the two phases, ra- 
diative cooling of the hot X-ray gas to form the cooler 
gas, and mixing of the hot and cool gas. Directly to the 
north and south of the galaxy center, there are holes in 
the X-ray emission. These holes are roughly coincident 
with the lobes of the radio source, suggesting that there 
may be two radio bubbles being blown in the hot gas. Fi- 
nally we sec that NGC 1603, a galaxy east of NGC 1600, 
has a tail of soft X-ray emission to its west. Calculations 



indicate that this is likely to be a ram-pressure stripped 
tail of ISM from NGC 1603, removed as a result of mo- 
tions of this galaxy through the surrounding group gas. 
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Discrete X-ray Sources in NGC 1600 
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TABLE 1 — Continued 
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R.A. 


Dec. 


d 




Count Rate 














No. 




Name 




(h m s) 


(° ' ") 


(") 


(") 


(10-'' 


S/N 


Lx 


H21° 


H31° 


H32° 


Notes 


(1) 




(2) 






(3) 


(4) 


(5) 


(6) 
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Note. — The units for Lx are 10^^ erg in the 0.3-10 keV band. 

^ Sources near the center may be confused with nearby sources, making their positions, fluxes, and extents uncertain. 
^ Source is noticeably more extended than PSF. 

^ Source is at the edge of the S3 detector, and flux is uncertain due to large exposure correction. 

^ Possible optical counterpart. 

® Source may be variable. 

^ Source is part of analysis sample. 



TABLE 2 
X-RAY Spectral Fits of NGC 1600 
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31+0-46 


1.11 


778 


76.3/65^1.17 



Units are 10"^^ erg cm-^ s~^ in 0.3-10 keV band. 

The adopted best-fit model for this emission. 

The energy range for this spectrum excludes 1.9 - 2.1 keV. 

The abundances for this spectrum were tied together. 



